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Abstract. Understanding the role of microenvironment in cancer growth and metastasis is a key issue for cancer re¬ 
search. Here, we study the effect of osmotic pressure on the functional properties of primary and metastatic melanoma 
cell lines. In particular, we experimentally quantify individual cell motility and transmigration capability. We then per¬ 
form a circular scratch assay to study how a cancer cell front invades an empty space. Our results show that primary 
melanoma cells are sensitive to a low osmotic pressure, while metastatic cells are less. To better understand the exper¬ 
imental results, we introduce and study a continuous model for the dynamics of a cell layer and a stochastic discrete 
model for cell proliferation and diffusion. The two models capture essential features of the experimental results and 
allow to make predictions for a wide range of experimentally measurable parameters. 
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1 Introduction 

The microenvironment is a key factor in tumour development 
and progression. Its influence on the behaviour of cancer cells 
is to a great extent mediated by the composition, structure, and 
dimensionality of the extracellular matrix, the polymeric scaf¬ 
fold that surrounds cells within tissues. Research shows that 
the mechanical properties of tumour microenvironment can fa¬ 
cilitate or oppose tumour growth and dynamics, although this 
effect is poorly understood. Mechanical stresses such as those 
experienced by cancer cells during the expansion of the tumour 
against the stromal tissue have been shown to release and ac¬ 
tivate growth factors involved in the progression of cancer (T). 
Moreover, the stiffness of the matrix surrounding a tumour de¬ 
termines how cancer cells polarize, adhere, contract, and mi¬ 
grate, thus regulating their invasiveness d . Another possibil¬ 
ity is that mechanical stresses directly regulate the growth and 
death rates of cancer cells. Recently, Montel et al. have ex¬ 
plored this possibility, investigating the effect of a constant 
stress applied on cellular spheroids over long time scales by 
inducing osmotic pressure by a solution of dextran, a biocom¬ 
patible polymer which is neutral and can not be metabolised 
by mammalian cells 120. Using a similar method, some of us 
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recently reported that a constant low osmotic pressure (lkPa) 
affects more strongly the proliferative capability of primary hu¬ 
man melanoma cells (IgR39) in comparison to the correspond¬ 
ing metastatic ones (IgR37) n . Furthermore, a computer sim¬ 
ulation analysis of the growth of melanocytic nevi inside epi¬ 
dermal or dermal tissues shows that the shape of the nevi is cor¬ 
related with the elastic properties of the surrounding tissue l4l. 

Several studies in the literature reported important changes 
in cellular functioning due to osmotic pressure 00, but the 
stresses involved (in the MPa range) were orders of magnitude 
larger than those (in the kPa range) studied in Refs. GIEIH. 
It is interesting to notice that compressive stresses of slightly 
less than 1 kPa applied through a piston were recently found 
to induce a metastatic phenotype in cancer cells m. In a re¬ 
cent paper, Simonsen et al. showed that interstitial fluid pres¬ 
sure (IFP) 0| was associated with high geometric resistance 
to blood flow caused by tumour-line specific vascular abnor¬ 
malities in xenografted tumours from two human melanoma 
lines with different angiogenic profiles j9j. In another recent 
paper, Wu et al. investigated how nonlinear interactions among 
the vascular and lymphatic networks and proliferating tumour 
cells might influence IFP transport of oxygen, and tumour pro¬ 
gression fna . They also investigated the possible consequences 
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of tumour-associated pathologies such as elevated vascular hy¬ 
draulic conductivities and decreased osmotic pressure differ¬ 
ences. All these parameters might affect microenvironmental 
transport barriers, and the tumour invasive and metastatic po¬ 
tential, opening interesting new therapeutic perspectives. The 
role of IFP for drug delivery was also analyzed computationally 
in Ref. CD . In general, understanding the influence of mechan¬ 
ical stress on cancer growth could shed new light on tumour 
development and progression. Studies present in the literature 
also use hypertonic conditions, but the effects are always ob¬ 
served at high osmotic pressure [ 12, 13]. 

In this paper, we show that low osmotic pressure (lkPa) 
besides proliferation, which was the focus of previous works 
ill also modulates key functional biological aspects of 
cancer cells such as their motility and transmigration, possibly 
contributing to the acquisition of a more aggressive phenotype. 
The choice of an osmotic pressure of lkPa was suggested by 
previous experiments (4) showing that this is the smallest value 
that has an appreciable effect on the growth of melanoma cells. 
Melanoma is an aggressive, radio- and chemo-resistant tumour 
which becomes impossible to cure when metastasised. There¬ 
fore, a pressing clinical problem in the treatment of melanoma 
is to understand how to stop or prevent the capability of the 
tumour to give rise to metastasis. Using an interdisciplinary 
approach combining cellular biology and theoretical physics, 
we show that a low osmotic pressure (lkPa) leads to significant 
changes in cell F-actin organization compromising the capabil¬ 
ity of the cells to move and transmigrate. Interestingly, osmotic 
pressure is more effective on the primary human melanoma cell 
line with respect to the metastatic one belonging to the same 
patient. We perform quantitative analysis on the experimental 
data and extract parameters for cell motilities that we can then 
use in a continuum theory and a discrete model for cancer cell 
front dynamics. Overall, our results suggest that osmotic pres¬ 
sure could contribute to the selection of a subpopulation in the 
primary human melanoma cells contributing to the acquisition 
of a more aggressive phenotype. 


2 Experimental results 

2.1 Effect of osmotic pressure on F-actin 
organization 



Fig. 1 Effect of osmotic pressure on F-actin. Subconfluent cells 
were plated on 33cm 2 Petri dish untreated or pre-coated with type I 
collagen. After 6 days the cells were fixed with 3.7% paraformaldeide 
and incubated with Red-Falloidin for 30min at room temperature. 

The nuclei were counterstained with DAPI and the slides were 
mounted with Pro-long anti fade reagent (Invitrogen). The images 
were acquired using a Leica TCS NT confocal microscope. Panel a-d 
represent IgR39 cells and panel e-f IgR37 cells. Scale bar is 50/im 


2.2 Effect of osmotic pressure on transmigration 
capability of lgR39 and lgR37 human melanoma cells 

An important function needed by cancer cells in order to metas- 
tasise is the ability to bypass an obstacle such as an endothelial 
wall. To simulate this process in vitro, we use a transwell as¬ 
say counting the number of cells that are able to pass through 
8/im wide pore after 6 hours in presence of a chemoattractant 
(10%FBS). Without chemoattractant the cells barely overcome 
the obstacles, although metastatic cells are intrinsically more 
capable to do this (Table [T}. This effect should, however, be 
considered just as background noise. Figure [2] shows a sig¬ 
nificant decrease of transmigration capability for cells that are 
treated six days with a lkPa osmotic pressure. We expose cells 
for some days to osmotic pressure in order to avoid measuring 
transient effects. The decrease is by a factor 6 for IgR39 and by 
a factor 2 for IgR37 with respect to the untreated ones (Table 
[l}. These results show that osmotic pressure reduces the trans¬ 
migration capability of melanoma cells but the effect is much 
stronger in primary IgR39 cells than in metastatic IgR37 cells. 


The organization of F-actin in the two human melanoma cell 
lines is shown in Figure 1 under normal growth condition or 
lkPa osmotic pressure without or with collagen as physiologi¬ 
cal substrate. In both cell lines, lkPa induces a rearrangement 
in F-actin organization (Fig. [TJ. In particular in IgR39, lKPa 
osmotic pressure induces the appearance of filopodia and stress 
fibers (Fig.[l^-b) Similar changes have been shown plating the 
cells on a collagen pre-coated Petri dish (Fig.JlJ-d). Indeed col¬ 
lagen should create a resistance against cell deformations that 
is similar to a confining pressure. The metastatic cell line IgR37 
shows, however, is only slightly affected by osmotic pressure 
displaying a more elongated shape (Fig[]J-f). 


2.3 Effect of osmotic pressure on circular scratch 
assay of lgR39 human melanoma cells. 

We study the biological effect of lkPa osmotic pressure on 
the capability of the cells to cover a free space using a circu¬ 
lar scratch assay. As in Fig. [3] the circular shape of the disk 
is maintained during up to 96 hours, suggesting that the cells 
move isotropically into the free space. Moreover, in Table [2] is 
shown that osmotic pressure induces a little effect on the migra¬ 
tion capability to cover the free area (around 4%). By observing 
the cells at higher magnification, we notice that the structure of 
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IgR39 OkPa 

IgR37 OkPa 

IgR39 lkPa 

IgR37 lkPa 

without FBS 
+FBS 

103 ±14 
1704 ±388 

602 ±74 
1925 ±173 

68 ±11 

269 ±21 

387 ±57 
825 ±125 


Table 1 Result of the transwell assay for primary and metastatic human melanoma cells with and without osmotic pressure. The results are 
the mean number of cells crossing pores averaged over three independent experiments. Error bars correspond to the standard error. 


Time [hours] 

Pressure [kPa] 

Empty area [mm 2 ] 

% of area covered 

0 

0 

32.1 ±0.2 

0±1 

24 

0 

21 ±2 

35 ±5 

24 

1 

23 ±5 

30± 15 

48 

0 

16 ± 1 

51 ±4 

48 

1 

21 ±4 

35 ±10 

72 

0 

17 ± 1 

47 ±3 

72 

1 

14±5 

57 ±16 

96 

0 

9±3 

70 ±9 

96 

1 

12 ± 1 

61 ±2 


Table 2 The table shows the free area and the percentage of covered area obtained from Fig. [5] Data are averaged over two independent 
experiments and reported with the associated standard deviation. 


- FBS 


a) 

b) 


OkPa 

lkPa 



10% FBS 

C) 

d) 

. : ■ ■ 

OkPa 

lkPa 




Fig. 2 Transwell Assay. 70000 cells are plated to each transwell 
insert (8/im pores) into the upper chamber. In the lower chamber 
600^1 medium without FBS (-FBS) or when 10% FBS is added. 
After 6 hours the cells are fixed and stained in a 20% methanol/0.1% 
crystal violet solution for three minutes at room temperature, 
followed by washing in deionized water to remove redundant 
staining. Non-migrated cells remaining at the upper side of the 
membranes, are carefully removed with cotton swabs and inserts are 
dried in darkness overnight. Pictures of the entire Petri dish are made 
with a Leica MZFL11 microscope mounted with a camera Leica 
DFC 32. The blue cells are counted using the magnified images 
(25X) and then calculated for the whole surface (32cm 2 ). 

the front is only slightly affected by osmotic pressur^3]>c.We 
also quantify the front roughness and find no significant differ¬ 
ence for different osmotic pressures or times. Colony formation 
assays reported in Ref. n for the same cells, show that the rate 
of division per day decreased from k^ w = 0.60 to = 0.51 
under the effect of lkPa of osmotic pressure. Here we investi¬ 
gate if the proliferative capacity of the cells are different at the 


Fig. 3 Circular scratch assay. Cells were plated on 33cm 2 Petri dish 
with at the center a disk (32mm 2 ) at confluence.The day after, the 
disk was removed and cells were allowed to keep migrating and 
dividing for different times, before they were fixed 

with 3.7% paraformaldeide and stained with hematoxillin/eosin 
solution at different times. The same experiment was carried out on 
cells treated with lkPa osmotic pressure. Panel a shows a typical 
experiment at magnification 8X; Panel b-c show details from the 
same images at magnification 40X. The images were acquired with 
LEICA MZFLIII mounted with a camera LEICA DFC320. 


border of the disk or at the edge of the Petri dish. As shown 
in Fig. [5j no significant differences of proliferation have been 
shown using Ki67, a known proliferative marker, between the 
centre of the dish and the edge (Fig. [5] shows the expression at 
72 hours, similar results was obtained at 24 and 48 hours). We 
did not carry out this kind of experiment in IgR37 since these 
cells do not grow as a monolayer and it is impossible to get a 
homogeneous confluent layer. 
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Fig. 4 Estimate of the front roughness, a) Plates resulting from the 
circular scratch assay are digitized and b) thresholded using the same 
exposure and the same threshold for all the images, c) The cell front 
is obtained by a percolation cluster algorithm and d) the roughness of 
the front is computed for each experimental condition (results of two 
independent experiments). For comparison, we also report the 
roughness obtained from the simulations of the computational model. 


2.4 Osmotic pressure affects cell motility 

To investigate the effect of osmotic pressure on cell motility, 
we reconstruct cell trajectories from confocal images (see [6^). 
We use IgR39 cells since these cells are the ones that are most 
affected by osmotic pressure (4) . To quantify cell motility we 
compute the typical width of the trajectories W = ((R(t) — 
( R )) 2 ), where the average is taken over all time frames and 
4 different cells for each experimental condition (see methods 
section). Fig. shows that osmotic pressure reduces the av¬ 
erage excursions of the cells but due to the large fluctuations 
the result is not statistically significant.To better differentiate 
between the two cases, we compute the time dependence of 
the mean square displacement ((R(t + t f ) — R(t')) 2 ) which in¬ 
creases roughly linearly in time (Fig. [6]c) as would be expected 
for a Brownian motion. The slope of the linear part of the curve 
can be used to estimate the diffusion coefficient resulting in 
D = 0.016(/im) 2 /s at OkPa and D = 0.004(jUm) 2 /s at lkPa. 
To confirm the Brownian nature of the process, we compute 
the power spectrum S(co) of the trajectories, showing that it 
decays approximately as 1 /ft) 2 as expected for a random walk. 


3 A continuous model for advancing cell 
front 

Since the main objective of our work is to study the effect 
of mechanical stress on an layer made of cancerous cells, we 
propose a physical model to explain and quantify the obser¬ 
vations. As shown in Table [2j we get a strong decrease of the 
area of the hole, about 50%, the first 2 days, then a slowest de¬ 
crease as shown in Fig. 0. The effect is a little less acute if an 
over-pressure is added via addition of dextran. For an advanc¬ 
ing layer under chemotaxis, a continuous model has already 



Fig. 5 Ki67 Immunohystochemistry IgR39 cells were plated on 
33cm 2 Petri dish with at the center a disk (32mm 2 ) at confluence. 72 
hours after the disk was removed, cells were fixed with 3.7% 
paraformaldeide and stained overnight with anti-mouse Ki67 (1:100, 
DAKO). After a brief incubation for 30min with biotinylated 
anti-mouse secondary antibody(Dako, 1:100), the cells were 
incubated with streptavidin (Dako, 1:350) conjugated to peroxidase 
for 30min. Color was routinely developed using DAB peroxidase 
substrate kit (Vector) up to 10 min and coverslipped with a 
permanent mounting medium. 


OkPA 


1 kPA 


t 






Fig. 6 Cell motility. Subconfluent cells plated in a 35cm 2 Petri dish 
untreated or treated for 6 days with lkPa of osmotic pressure are 
transfected with H2B-GFP plasmide. Pictures are made on a confocal 
microscope (Leica TCS SP5 AOBS with resonant scanner, equipped 
with a 20X, 0.5 NA, dry objective) every 4 minutes for 14 hours, a) 
Typical trajectories for IgR39 cells untreated and under osmotic 
pressure, b) The width of the trajectories is reduced under osmotic 
pressure but the result is not statistically significant due to the large 
error bars (n=4 for cells for each cases where considered), c) The 
average displacement grows roughly linearly with a typical slope that 
depends on the osmotic pressure, d) The Brownian nature of the 
motion is confirmed by the power spectrum of the trajectories that 
decays as ft) -2 . 


been proposed mma , but these works are not concerned by a 
proliferative layer submitted eventually to an overpressure due 
to an osmotic solution. In addition, proliferation is not often 
considered in the advance of a monolayer but here our experi¬ 
ment clearly shows that it cannot be discarded. Let us consider 
a layer of constant thickness in a domain Q with a border 8Q. 
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The cell displacement in £2 is mostly driven by the cell pro¬ 
liferation, but at the border the cells may be more proliferative 
eventually and even detach and move freely. Due to prolifera¬ 
tion, the mass conservation equation is then transformed into: 

^ +V(poV) = -k v (P-P h ), (1) 


where P/ ? is the homeostatic pressure for which cell prolifera¬ 
tion just compensates apoptosis, po is the constant mass den¬ 
sity, V is the velocity at an arbitrary point of Q. This sim¬ 
ple idea about regulation in living systems is due to Claude 
Bernard, a french physiologist, in 1865 on. Mathematically, 
such a continuous model is valid at a scale larger than the cell 
size and the velocity is an averaged quantity on a sample of 
intermediate scale between the single cell and the full layer. 
Adding dextran to the solution increases the pressure P, slightly 
inhibiting the proliferation rate. For a layer crawling on a sub¬ 
strate in strong adhesion, the second Newton’s law of dynamic 
simplifies into a Darcy law ([17, 14]) 


y = —vp, 

M 


( 2 ) 


where the mobility coefficient M is related to the friction coef¬ 
ficient. Combining both equations gives: 

AP- a 2 P = 0 with P = P-P h (3) 

and a 2 = (k v MRf)/po, with Rj the initial radius of the hole of 
order the millimeter, chosen as our length unit in the follow¬ 
ing. At the free interface, boundary conditions include the me¬ 
chanical equilibrium and continuity of normal velocities. The 
pressure at the layer border Pj is equal to the pressure of the 
solution (if capillary effect is neglected). It is smaller at the in¬ 
terface S Q than in the layer bulk Q for a proliferative layer. 
The interface velocity S£2 t is given by the velocity of the cells 
at the interface which may be different from the cells in the 
bulk Q . Indeed, these peripheral cells have more room both for 
displacement and proliferation. Then the normal front velocity 
8Q t is given by 


N8Q t = -—N-VP + N-v s (4) 

M 

where may represent a specific proliferation rate of the bor¬ 
der cells. Other boundary conditions may exist due to the ge¬ 
ometry of the experiment, that is why we consider now 3 dif¬ 
ferent growth geometries. 


3.1 Growing circular disc 


We imagine a disc of a cell monolayer and assume that the 
growth process preserves the circular geometry. Eq. 0 for the 
pressure P = P — P^ is then transformed into an ordinary differ¬ 
ential equation of second order: 


1 d 
RdR 


f dP\ 9 

(R ^ =ae 


(5) 


The regular solution at the origin is the modified Bessel func¬ 
tion of zero order 7 q so 


P(R) = Pi 


IpjaR) 

Io{aRi(t )) 


(6) 


and Pi = Pi— Ph is negative, being fixed by the nutrient medium, 
whatever the hole radius Ri(t). From Darcy’s law (Eq. ([2])), we 
derive the growing velocity at the interface: 


N-8Q t =Ri(t ) 


p a h(aRj(t)) 

‘MI 0 (aRi(t )) 


(7) 


This equation means that the interface velocity has the same ve¬ 
locity as the moving cells of the interface. As for the Saffman- 
Taylor viscous fingering in radial geometry m, our set of 
equations is time independent but the time dependence reap¬ 
pears via the condition of continuity of velocities given by Eq. 0 
It is an implicit equation for the time variation of the radius 
which can be solved numerically. However one can make an 
asymptotic analysis of Eq. 0 and gets a constant growth rate at 
long times : 

and a radius growing linearly in time. One can notice that a con¬ 
stant proliferation rate into Eq.Q will lead to an exponential 
time behavior for the radius so to follow the growth expansion 
of the layer at long times may be a good test for the validation 
of the model. 


3.2 A hole in an infinite layer 


In this case, we can accept the irregular Bessel function at the 
center Kq as the solution of Eq.Q, but not /o, which grows 
exponentially. The solution reads 


P(R) = Pi 


Ko(ccR) 

Ko(aRi{t))' 


( 8 ) 


The pressure P vanishes at infinity where there is no cell pro¬ 
liferation. The velocity given by Darcy’s law is then: 


O' KijaR) 
M l K 0 (aRi(t)) 


(9) 


being negative as Pi and vanishing as R -A oo as expected. Here 
again the radius of the hole is given numerically if we take 
Eq.Q on the interface. In addition, the radius velocity diverges 
when the total cover of the hole is reached, and the continu¬ 
ous model stops to be valid for Ri(t) -A 0. Indeed, Eq.([9]) is 
a nonlinear first order differential equation for the hole radius 
that we cannot integrate explicitly. However, for a large enough 
hole, assuming a of order 1 , the decrease of the radius is linear 
at initial time, but the radius will finally vanish like a square 
root, (t — t c ) 1 / 2 , near the complete closure arising at time t c . 
Before presenting more elaborated model for our experiment, 
we would like to comment about the stability analysis of the 
border and the assumption of circular symmetry for our set of 
equations. A careful observation of our experiments presented 
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in Fig. 0 shows that the inner border is somewhat noisy but 
does not present really instabilities similar to patterns observed 
in viscous flow ED or bacteria colonies Da. It turns out that, 
even for a planar front, this model is stable against transverse 
perturbations, even in absence of capillary effect. Although as¬ 
tonishing, it is in agreement with observation. In order to con¬ 
vince the reader, we simply present the stability treatment for 
this geometry in the Appendix section and analytically prove 
that our model predicts that the migration induced by growth 
is stable in the circular geometry. So the border irregularities 
observed in Fig. ^ are not fingering but can be interpreted as 
the stochastic behavior of tumor cells. 



Fig. 7 Front propagation. Relative area of the hole as a function of 
time for experiments (see Table[2]), theory and simulations. 
Parameters for the theory are P[ = — 1 for OkPa and P[ = —0.98 for 
lkPa, (Xi = 0.71, a e = 0.53, 8 = 0.3 (with initial radius 1), z = 1. 
These numerical values correspond to the best estimated fitting for 
our experimental data. Parameters for the simulations are = 32 
division/day while keeping &diff Adiv = 0.9 or 1.7, for OkPa and lkPa, 
respectively. 


3.3 The circular scratch assay 

Here we focus on the geometry of the circular scratch assay: 
it is made of a planar annulus of cells limited by the rigid cir¬ 
cular wall border of the Petri dish at R = R e and , at R = Ri, 
is free to move into the empty space. The interior interface at 
R = Ri is rather diffuse but an average circular border can be 
defined. In addition a thin circular ring (of thickness 8) with 
smaller cell density can be observed appearing less contrasted 
than the inner core, located at R = Ri + 8. We will adapt the 
previous model to this inhomogeneous ring which has a den¬ 
sity pi for R between Rj(t) and Ri(t) + 8 while the density in 
the bulk (Ri(t) + 8 < R < R e ) is called p e . It is important to 
remember that p is indeed a density per unit surface, the volu¬ 
metric density being the same for all cells. As a consequence, if 
we maintain the hypothesis of a monolayer, the surface density 
may be larger at the inner border where the cells are more free 


to advance than in the annulus where they are more confined. 
It explains the optical difference in Fig. 0- We take constant 
both p values. Moreover, in vitro, it has been shown that in 
confluent layers, live cell extrusion is induced to recover the 
homeostatic state in case of over-crowding and it is clear that 
the border allows more easily such extrusion l20l . So we pre¬ 
sume that pi < p e . Because the layer geometry excludes R = 0 
and R -A oo, possible solutions for Eq. 0 are superpositions of 
Iq and Kq and we get 


A (R) 


Pj_ r k(OjR) Ko(aiR) \ 

1 + b \ I(){a t R t (t)) + KoicCjRiit)) J 


( 10 ) 


with af = k v M/pj. The velocity of the cells at the outer border 
being : 


Oj Pi f hjOiR) K t (a,R) -j 

M\+b\l {) (a,R,(t)) KoiCkRiit)) f 

Pi keeps the same definition as before : the pressure in the so¬ 
lution minus the homeostatic pressure. For Ri(t) + 8 < R < R e 
we have: 


P e {R)=P e { 


kjcCeR) 
I\ ( K dc e R e ) 


KojOCeR ) \ ^ 

K\(a e R e ) J 


( 12 ) 


T being an evaporation rate of cells as mentioned by 120) and 


otv p f h(tt e R) K\{a e R) \ 
M e \h{a e R e ) K\ (a e R e ) I 


with oig = (k v MRj)/p e . In Eqs.(12 131, we have applied the 
cancellation of cell velocity at the exterior border R e . We have 
2 unknowns: the coefficient b and P e which are fixed by the 
mechanical equilibrium P l (Ri(t ) -j- 8) = P e (Ri(t ) + 8) and the 
continuity of the flux: PidP 1 /dR = p e dP e /dR for R = Ri{t) + 
8. The principle of such analysis is simple but requires also 
a numerical solution that we present in Fig.Q in compari¬ 
son with experimental data. We have four unknown parame¬ 
ters a z , a e , T and Pi. Varying parameters, the comparison with 
the experiment is mostly controlled by the quantity aPi/M = 
\Jk v /(piM)Ri which is the inverse of a typical time. In the ex¬ 
periment the time scale is the day. Taking into account the ex¬ 
perimental results shown in Fig. ([ 7 ]), the relative area variation 
of the hole after 2 hours is 0.1 for an osmotic pressure varia¬ 
tion of 1 kPa. It corresponds to a variation of the front veloc¬ 
ity in international unit of 710 -9 m/s (evaluated from the ex¬ 
perimental data of Fig. ([ 7 ])) which gives the unknown value of 
M of 10 14 Pas/m 2 , using Eq(T3]). Once we have this value we 
can derived easily the proliferation rate k v /po of order 510 -9 
SI. The friction coefficient M for an advancing epithelium has 
been evaluated in 03 for MDCK cells giving a value an order 
of magnitude larger 10 15 . Of course, the value of M depends on 
the cells involved in the experimental set-up, and are not uni¬ 
versal. However this estimation is consistent with usual order 
of magnitude for biophysical values. Taking again the order of 
magnitude of the velocity as 10“ 8 m/s, knowing that forces one 
needs or can apply to a cell are typically of order pico-Newton 
so 10“ 12 N gives a friction coefficient of 10 14 Ns/m 4 (being the 
ratio of the force value by cell volume and by cell velocity). Fi¬ 
nally sensitive variation by osmotic pressure being 0.02 in our 
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model which corresponds to 1 kPa, gives a naive estimation of 
50 kPa for the departure from the homeostatic pressure. A bet¬ 
ter calculation can be done from the model using Eq.fTT]) and 
gives 25 kPa. 


4 Computational model of cell duplication 
and movement 

To understand how the morphology of the front in the circu¬ 
lar scratch assay depends on the interplay between cell divi¬ 
sion and mobility, we simulate a simple lattice model includ¬ 
ing mechanisms for cell division and motility, with the rela¬ 
tive rates of the two as the tunable parameter. The model is de¬ 
fined on a square lattice of linear size L = 700, in which all the 
sites are occupied by cells except those inside a circle of radius 
r = 320. Cells that have at least one unoccupied nearest neigh¬ 
bour site can either divide with rate or diffuse with rate 
ha iff. The kinetics of the model follows the Gillespie algorithm 
ED . At each time step we count how many cells have 1,2, 3, or 
4 neighbours (i.e. n\,ri 2 , n 3 and ^ 4 ). Then the total rate is taken 
to be Rt 0 t = (n\ + 0.75 *^3 +0.5 *^2 + 0.25 *n\) * (&div + &diff)> 
where the reaction rate is taken to be proportional to the num¬ 
ber of neighbours a cell has, with the bare rates given by those 
of a cell with all its four nearest neighbours vacant. At each 
time step, the reaction to take place is chosen randomly, with a 
probability (rate of the reaction)/7? tot . There are 8 possible reac¬ 
tions with different rates, i.e. 4 possible division reactions and 
4 possible diffusion reactions, each corresponding to a different 
number of nearest neighbours a cell has (1,2,3 or 4). First, the 
reaction to take place is chosen, and then a random cell with the 
number of nearest neighbours corresponding to the selected re¬ 
action is found, and the reaction is carried out. The correspond¬ 
ing time step is generated from an exponential distribution with 
R to t as the rate parameter. Then the simulation proceeds to the 
next time step, and the process continues until cells invade the 
center of the hole. Results are averaged over 10 realizations 
for each case. In Fig. [ 8 ] we report typical configurations of the 
circular scratch assay obtained by numerical simulations for 
different values of fcdiffAdiv- When this ratio is large the in¬ 
terface is very diffuse, with several cells moving ahead of the 
front. As the ratio decreases, the front is more sharp, but still 
remains rough due to the randomness in the division process. 
In order to compare our simulations with experiments we first 
estimate the ratio fcdiffAdiv from other experiments. In partic¬ 
ular, we can estimate k^ff ~ D/a 2 , where a is the typical cell 
diameter. From Fig. [I] we obtain that a ~ 20/im at OkPa and 
a ~ 50/im at lkPa. Using our measurements for D and 
0 , we obtain kdiff Adiv — 0.9 at OkPa and kdiff Adiv — 1-7 at 
lkPa. If we use these value of k^ v and &diff, the simulated front 
advances too slowly with respect to the experiments. This is 
expected since in the circular scratch assay cells are induced to 
rapidly fill the hole. A good agreement between model and ex¬ 
periments is obtained setting &div = 32 division/day while keep¬ 
ing &diff Adiv = 0.9 or 1.7, for OkPa and lkPa, respectively. The 
result is reported in Fig. [ 7 ] for the front evolution and in Fig. 
[4] for the roughness. The parameters employed in the model to 
obtain the best fit to the experiments look, however, unrealistic 
when compared with estimates obtained from experiments on 


individual cells. One possibility could be the presence of ac¬ 
tive forces on the cells that assists the closure of the hole. In 
the model, we have also quantified the front velocity and the 
roughness that are reported in Fig. [9] The velocity is found to 
decrease when diffusion predominates over cell division, sug¬ 
gesting that the process of ring closure is dominated by cell di¬ 
vision (Fig. [9^). Fig. [ 9]3 shows instead that the front roughness 
increases when diffusion becomes more important. The effect 
of osmotic pressure is, however, very small and is hard to detect 
by looking at the roughness, in agreement with experiments. 



Fig. 8 Simulated ring assay. Illustrative examples of the rings 
obtained for different ratios of the rates of diffusion and division. 



Fig. 9 Front velocity and roughness from numerical simulations. 

a) The front velocity as a function of the ratio of the rates of diffusion 
and division, b) The front roughness as a function of the ratio of the 
rates of diffusion and division. 
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5 Discussion 

Melanoma is one of the most aggressive tumours and practi¬ 
cally impossible to cure once it becomes metastatic. Hence, 
one of the main goals of research is to understand the processes 
involved in metastasis. The classical approach is to consider 
the biological properties of the cell trying to identify biologi¬ 
cal markers. Here, we investigated this problem from a differ¬ 
ent point of view, focusing on the possible role of mechani¬ 
cal stress in the selection of a more aggressive cell subpopula¬ 
tion. Recent research has shown the involvement of mechani¬ 
cal stresses, such as the osmotic pressure, in the modulation of 
critical biological properties EE [3IU El. Therefore, we have 
investigated the effect of low osmotic pressure on two cell lines 
belonging to the same patient, IgR39 and IgR37 a primary and 
metastatic cell line, respectively. Our results show that a low 
osmotic pressure which gives little effect on cell proliferation 
on the primary human melanoma cell line si, induces signif¬ 
icant changes on the F-actin organization,with the appearance 
of filopodia and stress fibres. These changes are more evident 
in primary human melanoma cell lines than in metastatic ones. 
Earlier studies on endothelial cells have shown that an osmotic 
pressure around 1.6 kPa induces an elongation of the cells due 
to a rearrangement of the cytoskeleton l22ll23l . in agreement 
with our results. In human breast cancer, mechanical compres¬ 
sion has an effect on the contractility of the cells suggesting the 
presence of a mechano-regulation which is able to integrate ex¬ 
ternal physics stimulation into cytoskeletal changes and there¬ 
fore cell phenotypes 0- A critical question is whether F-actin 
rearrangement has a functional effect on biological properties 
of cancer cells such as their capability to move or to invade. 
Metastasis is indeed an integrated process composed by sev¬ 
eral steps. At least two key aspects play a critical role in this 
process: first the capability of the cell to move and then the 
capability of the cell to squeeze through the endothelium. To 
reproduce in vitro these two processes we have carried out the 
circular scratch assay which gives the possibility to study the 
movement of the cell into an open space and the transwell assay 
which investigates the capability of the cell to squeeze through 
a pore of a fixed size. Primary human melanoma cell are less 
able to move through fixed size pores using transwell assay 
in comparison to metastatic cells. Moreover, interestingly, low 
osmotic pressure leads to a reduced capability to transmigrate 
in both cell lines, however, the primary human melanoma cell 
line is more affected. The circular scratch assay confirms an 
inhibitory effect of osmotic pressure on the capability to move 
of the cells. This is also confirmed by quantitative measure¬ 
ments of cell motility from series of confocal images. In or¬ 
der to better understand the biophysical aspects of the circular 
scratch experiments, we have approached the problem theoreti¬ 
cally and numerically. We have derived a continuum theory for 
the growth of a single layer of cells with different boundary 
conditions such those used in the circular scratch assay. This 
allows to obtain quantitative predictions for the dynamics of 
the cell front. In particular, we calculate the dependence of the 
front velocity on a few experimentally measurable parameters 
such as the cell division rate and motility. The theoretical re¬ 
sults are in good agreement with the experiments. The main 
limitation of the theoretical model is that it does not include 
noise and fluctuations which are obviously present in our ex¬ 


periments. To deal with this issue in a simple way, we introduce 
a discrete lattice model of cell proliferation and diffusion which 
we simulate in similar conditions. The discrete model allows to 
reproduce the experimentally observed roughening of the cell 
front. This combination of theoretical and numerical models al¬ 
lows to understand better the role of motility and cell division 
under osmotic pressure and provide a useful tool to explore the 
effect of experimentally relevant parameter. Both cell prolifer¬ 
ation and motility rates are reduced by osmotic pressure, which 
also changes the shape of cells making them thinner and more 
elongated. The combination of this factor leads to an slight in¬ 
crease of &diff Adiv under osmotic pressure. All together our re¬ 
sults show how low osmotic pressure changes functional prop¬ 
erties of tumour cells but have less effect on more aggressive 
metastatic cell populations. These results might be in agree¬ 
ment with previous papers showing a role of osmotic pressure 
in the selection of more resistant subpopulation l24l fll 251. 


6 Material and Methods 

6.1 Dextran solution 

A master solution of dextran at 10% (w/v) was formed (Dextran 
from Leuconostoc spp, Fluka) and diluted to the desired con¬ 
centration with complete medium. Transformation from dex¬ 
tran concentration to osmotic pressure was performed accord¬ 
ing to the calibration curve measured in Ref. [26). 

6.2 Cell lines 

Human IGR39 and IGR37 cells were obtained from Deutsche 
Sammlung von Mikroorganismen und Zellkulturen GmbH and 
cultured as previously described G2 . IGR39 was derived from 
a primary amelanotic cutaneous tumour and IGR37 from an 
inguinal lymph node metastasis in the same patient. The cells 
were plated on pre-coated collagen type I plates according to 
the manufacturer’s instructions for further experiments when 
specified in the results section. 

6.3 Cell transfection 

Cells were transfected with l fig H2B-GFP plasmide (plasmide 
116980, Addgene) according to the manufactories instruction 
of SuperFect Transfection reagent by Qiagen. The cells were 
analyzed by time lapse using Leica TCS SP5 AOBS with reso¬ 
nant scanner, equipped with a 20X, 0.5 NA, dry objective (Le¬ 
ica Microsystems GmbH, Wetzlar Germany) with fluorescence 
under 5% C02 and 37C temperature. 

6.4 Immunofluorescence-Phalloidin 

Subconfluent cells grown on glass coverslip were fixed with 
3.7% parafolmaldeide in PBS for lOmin, permealized with 0.5% 
Triton X-100 in PBS for 5min at room temperature and. stained 

TM 

with 200/il of lOOnM Acti-stain 555 Phalloidin (TebuBio, 
Cat. # PHDH1) for 30min. DNA was counterstained for 30s 
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with 200/11 of lOOnM DAPI and the slides were mounted with 
Pro-long anti fade reagent (Invitrogen). The images were ac¬ 
quired using a Leica TCS NT confocal microscope.The same 
experiments were carried out using pre-treated collagen type I 
coverslips (Sigma). 


6.5 lmmunohystochemistry-Ki67 


Cells plated for the circular scratch assay were fixed after 72 
hours with 3.7%paraformaldeide for lOmin at room tempera¬ 
ture, incubated for lhr in 1%BSA, 10NGS, 0.3M glycin, 0.1%Tween 
in PBS and overnight with anti-Ki67 (clone MIB-1, Dako, 1:100 6.8 Cell tracking 
cod.) at 4C. All the sections were incubated for 30min with bi¬ 
otinylated anti-mouse secondary antibody in PBS (Dako, 1:100). 

This was followed by incubation with streptavidin (Dako, 1:350) 
coniugated to peroxidase in PBS for 30min, and by a brief rinse 
in PBS. Color was routinely developed using DAB peroxidase 
substrate kit (Vector) up to 10 min, rinsed in distilled water, and 
coverslipped with a permanent mounting medium. The images 
were acquired using a microscope Leica MZ FLIII mounted 
with a camera Leica DFC320. 


for three minutes at room temperature, followed by washing 
in deionized water to remove redundant staining [16]. Non- 
migrated cells remaining at the upper side of the membranes 
were carefully removed with cotton swabs and inserts were 
dried in darkness overnight. The following day at stained mem¬ 
branes were made pictures using a transmitted-light microscope 
(Leica MZFL11 microscope mounted with a camera Leica DFC 
32) at different magnification (8, 16 and 25X). The blue cells 
were counted using the magnified images (25X) and then cal¬ 
culated for the whole surface (32cm 2 ). 


Cell tracking has been done in ImageJ by first projecting the 
confocal stacks (maximum intensity) and then finding the max¬ 
ima for different frames. In each frame i = 1, ...,V, we iden¬ 
tify the cell center of mass (X;,U) and then compute the width 
of the trajectory as W = £,((X; - X ) 1 + (Y, - Y) 2 )/N, where 
X = Y.,Xi/N and Y = £, Yj/N. The mean-square displacement 
is calculated as 


C(i) 


6.6 Circular scratch assay 


N/2 

= 1 2 m+j- 

j =i 


-Xj ) 2 + (Y i+ j- 


Yj) 2 )/N, 


(14) 


Cells were plated on 33cm 2 Petri dish with a disk at the cen¬ 
ter (Sigma cod. z370789, diameter: 6.4mm, area 32mm 2 ) to 
reach 70-80% confluence as a monolayer. When they were 
at confluence, the disk was removed. At different times after 
removal (from 4 to 96 hours), the cells were fixed with 3.7% 
paraformaldeide a different times (from 4 to 96hrs) and stained 
with hematoxillin-eosin. The images were acquired using a Le¬ 
ica MZFL11 microscope mounted with a Camera Leica DFC 
32 at different magnification (8, 16 and 25X). To quantify the 
free area and the correspondent percentage of covered area, we 
used Gimp (GNU Image Manipulation Program) to threshold 
the images. The front profile is obtained by using a simple clus¬ 
ter algorithm, as illustrated in Fig. [4] We then compute the area 
enclosed by the front and average the result over two indepen¬ 
dent realizations of the experiment under the same experimen¬ 
tal conditions. We also use the same algorithm to obtain the 
front roughness < 7 , defined as the standard deviation of the front 
radii. 


and the power spectrum by taking the discrete Fourier trans¬ 
form of ( Xi,Yi ) and then squaring the amplitudes. 


6.9 Statistical analysis 

Statistical significance was evaluated using the Kolomogorov- 
Smirnov non parametric test (see http://www.physics.csbsju.edu/stats/KS- 
test.html). 


6.10 Appendix 

Let us consider the stability analysis of the model of section 


3.2| For that we expand the pressure field and the inner border 
P(R,d,t) = P{R) + ep{R)cos(mQ)e nmt 


as 


Ri(0,t) = Ri(t) + £cos(m0)e‘ 


,£2 m t 


(15) 

(16) 


6.7 Transwell migration assay 

Migration experiments were conducted using a conventional 
24-well Transwell system (6.5 mm TranswellH (# 3422), Com- 
ingH, NY, USA) with each well separated by a microporous 
polycarbonate membrane (10 gm thickness, 8 gm pores) into 
an upper (insert) and a lower chamber (well) according to [16]. 
After 24 hours of serum deprivation, cells were detached, counted 
and resuspended in growth media without FBS to obtain equal 
cell densities (296000 cells/cm 2 ). A volume of 250 /il contain¬ 
ing 70000 cells was plated to each insert and 600 p 1 medium 
was added to the wells. As chemotactic factors, a medium con¬ 
taining 10% FBS was used. After 6 hours the cells were fixed 
and stained in a 20% methanol/0.1% crystal violet solution 


where P(R), the pressure calculated in the circular geom etry 
and the inner hole radius Rf ( t ) are given respectively by Eq.( 8 9) 


Hereafter we simplify the notations into R[. Replacing P by P 
into Eq.([3]) gives for p(R) the following equation: 


Rrp” + Rp - (a 2 R 2 +m 2 )p = 0 


(17) 


whose regular solution at infinity is the modified Bessel func¬ 
tion of m order K m (aR)\ 


p{R) = po 


KmjocR) 
K m (aRj) 


(18) 


Because P = P t at /?/, the linear variation of the pressure van¬ 
ishes at R = Ri giving po +P'(Ri) = 0. So we deduce the value 
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of po = aPiKi(aRi)/Ko(aRi). The continuity equation modi¬ 
fied by the border distortion is then 

n m = -^(F’ + p') (19) 


After some elementary calculation involving the modified Bessel 
functions we get 


fl m 


a 2 V f | KijaRj) / 1 -m \ ^ 

Ml Ko(aRi) \ aRi K m {aRi) )) K ’ 


Taking the 2 extreme limiting values Ri -A 0 or Ri -a ©o, we find 
Q m always negative. Indeed for Ri going to zero, the second 
term in the bracket is dominant and negative (Pi < 0 ) 


£2 m ^ 


Pii i-^ 2 ) 
2 M 


( 21 ) 


which shows that Q m is always negative when the hole ra¬ 
dius decreases. This stability linear analysis proves that, for the 
circular geometry of the hole closing, induced by cell prolif¬ 
eration, our modeling predicts a stable process. Although the 
equations have a strong similarities with Laplacian growth HD 
this model predicts that there is no fingering as time goes on. 


7 Acknowledgments 

This work is supported in part by AAP Physique and Cancer 
2012, project DERMA (MBA). SZ acknowledges the visiting 
professor programm of UPMC, the European Research Coun¬ 
cil Advanced grant SIZEFFECTS and the Academy of Fin¬ 
land FiDiPro progam, project 13282993. CAMLP acknowl¬ 
edges the visiting professor programme of Aalto University 
and support from PRIN 2010. LL and MJA are supported by 
the Academy of Finland through an Academy Research Fel¬ 
lowship (LL., Project No. 268302) and through the Centres of 
Excellence Program (MJA and L. L., Project No. 251748). 


References 

1. M.J. Paszek, N. Zahir, K.R. Johnson, J.N. Lakins, G.I. 
Rozenberg, A. Gefen, C.A. Reinhart-King, S.S. Margulies, 
M. Dembo, D. Boettiger et al., Cancer Cell 8(3), 241 
(2005) 

2. F. Montel, M. Delarue, J. Elgeti, L. Malaquin, M. Basan, 
T. Risler, B. Cabane, D. Vignjevic, J. Prost, G. Cappello 
et al., Phys. Rev. Lett. 107, 188102 (2011) 

3. F. Montel, M. Delarue, J. Elgeti, D. Vignjevic, G. Cap¬ 
pello, J. Prost, New Journal of Physics 14(5), 055008 
( 2012 ) 

4. A. Taloni, A.A. Alemi, E. Ciusani, J.P. Sethna, S. Zapperi, 
C.A.M. La Porta, PLoS One 9(4), e94229 (2014) 

5. B. Racz, D. Reglodi, B. Fodor, B. Gasz, A. Lubics, 
F. Galiyas, Jr, E. Roth, B. Borsiczky, Bone 40(6), 1536 
(2007) 

6 . M.B. Nielsen, S.T. Christensen, E.K. Hoffmann, Am J 
Physiol Cell Physiol 294(4), Cl046 (2008) 


7. J.M. Tse, G. Cheng, J.A. Tyrrell, S.A. Wilcox-Adelman, 
Y. Boucher, R.K. Jain, L.L. Munn, Proc Natl Acad Sci U S 
A 109(3), 911 (2012) 

8. D. Fukumura, R.K. Jain, J. Cell. Biochem. 101, 937 (2007) 

9. T.G. Simonsen, J.V. Gaustad, M.N. Leinaas, E.K. Rofstad, 
PLoS One 7(6), e40006 (2012) 

10. M. Wu, H.B. Frieboes, S.R. McDougall, M.A.J. Chaplain, 
V. Cristini, J. Lowengrub, J Theor Biol 320, 131 (2013) 

11. M. Welter and H. Rieger PLoS ONE 8, e70395 (2013) 

12. O. Bounedjah, L. Hamon, P. Savarin, B. Desforges, P.A. 
Curmi, D. Pastre, J Biol Chem 287(4), 2446 (2012) 

13. Z. Ignatova, L.M. Gierasch, Proc Natl Acad Sci U S A 
103(36), 13357 (2006) 

14. M. Ben Amar, Eur Phys J E Soft Matter 36(6), 64 (2013) 

15. M. Ben Amar, M. Wu, J R Soc Interface 11(93), 20131038 
(2014) 

16. C. Bernard, Introduction a V etude de la medecine 
experimental (Paris, 1865) 

17. A.C. Callan-Jones, J.F. Joanny, J. Prost, Phys Rev Lett 
100(25), 258106 (2008) 

18. L. Patterson, Journ. of Fluid Mech. 113, 513 (1981) 

19. M.B. Amar, J. Phys.I France 2, 353 (1993) 

20. G.T. Eisenhoffer, P.D. Loftus, M. Yoshigi, H. Otsuna, C.B. 
Chien, P.A. Morcos, J. Rosenblatt, Nature 484(7395), 546 
( 2012 ) 

21. D.T. Gillespie, Journal of Computational Physics 22(4), 
403 (1976), I 

22. A.D. Acevedo, S.S. Bowser, M.E. Gerritsen, R. Bizios, J 
Cell Physiol 157(3), 603 (1993) 

23. S.A. Salwen, D.H. Szarowski, J.N. Turner, R. Bizios, Med 
Biol Eng Comput 36(4), 520 (1998) 

24. Y. Boucher, J.M. Kirkwood, D. Opacic, M. Desantis, R.K. 
Jain, Cancer Res 51(24), 6691 (1991) 

25. E.K. Rofstad, S.H. Tunheim, B. Mathiesen, B.A. Graff, 
E.F. Halspr, K. Nilsen, K. Galappathi, Cancer Res 62(3), 
661 (2002) 

26. C. Bonnet-Gonnet, L. Belloni, B. Cabane, Langmuir 
10(11), 4012 (1994), 

27. R. Taghizadeh, M. Noh, Y.H. Huh, E. Ciusani, L. Sigalotti, 
M. Maio, B. Arosio, M.R. Nicotra, P. Natali, J.L. Sherley 
et al., PLoS One 5(12), el5183 (2010) 







